Evaluating matrix elements relevant to some Lorenz violating operators 
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Carlson, Carone and Lebed have derived the Feynman rules for a consistent formulation of non- 
commutative QCD. The results they obtained were used to constrain the noncommutativity param- 
eter in Lorentz violating noncommutative field theories. However, their constraint depended upon 
an estimate of the matrix element of the quark level operator (jl — m) in a nucleon. In this paper 
we calculate the matrix element of (ft — m), using a variety of confinement potential models. Our 
results are within an order of magnitude agreement with the estimate made by Carlson et al. The 
constraints placed on the noncommutativity parameter were very strong, and are still quite severe 
even if weakened by an order of magnitude. 



I. INTRODUCTION 

In the recent literature, there have been considered a 
number of ways to modify the structure of space-time 
which can have experimental consequences. In one of 
the most popular scenarios, space-time is considered to 
become noncommutative at short distance scales, with 
space-time coordinates satisfying a commutation relation 
of the following form § |, 1 § |, g § 



ar,ar 



i9» v , 



where is a position four-vector promoted to an op- 
erator, and 9^ v is a set of c-numbers antisymmetric in 
their indexes. The most striking effects of space-time non 
commutativity of the form (|l|) are the Lorentz violating 
effects appearing in field theories, which is a consequence 
of 9° l and e l ^ k 9 %1 defining preferred directions in a given 
Lorentz frame. 

Jurco, Moller, Schraml, Schupp and Wess || have 
shown how to construct non-Abelian gauge theories in 
noncommutative spaces from a consistency relation. Us- 
ing the same approach Carlson, Carone and Lebed Q 
have derived the Feynman rules for consistent formu- 
lation of noncommutative QCD and they have com- 
puted the most dangerous, Lorentz-violating operator 
that is generated through radiative corrections. They 
have found that at the lowest order in perturbation the- 
ory, the formulation of noncommutative QCD that they 
have presented leads to Lorentz violating operators such 
as § 

O^qcr^q, Q^qo^pq and 0"" D^q<J up DPq. (2) 

In Q the phenomenological implications of the first of 
these operators were studied in detail. Noting that con- 
tributions from the space-space part of 9^ v make a llv ^ LV 
act like a a ■ B interaction with B directly related to 
9 lJ , a limit was placed on the scale of non commutativ- 
ity. One used the result of tests of Lorentz invariance in 



clock comparison experiments j8| , which suggest that ex- 
ternal a ■ B like interactions are bounded at the 10 _7 Hz 
level or few x 1CP 31 GeV. Carlson et al. || concluded 
that 



OA < 10 



-2!) 



(3) 



where 9 is a typical scale for elements of the matrix 9^. 

However, the effective Lorentz violating operator was 
obtained from a one loop correction to the quark propa- 
(1) gator, and the operator proportional to cr M „6' A " y also con- 



taincd a factor (jf — m). With B constant, the evaluation 
of a ■ B factors out from the evaluation of — m), and 
our discussion is focused on the later. 

In j| an ad hoc estimate was used for the matrix el- 
ement of the operator (ft — m), where m is the current 
quark mass, in getting limit in equation (||). The matrix 
element (ft — ni) was estimated to be about Mjv/3 « 300 
MeV, where Mjv is the nucleon mass. However, it has 
been argued that the expectation value of (ji— m) could 
be much less than this naive estimate H . 

The aim of this paper is to calculate the matrix element 
of the operator (ji — m), using variety of confinement po- 
tential models, so to evaluate the quality of the estimate 
made in 

The sample of potentials included four different con- 
fining potentials, two of them purely Lorentz scalar and 
two of them equal mixture of scalar and vector. The first 
scalar potential is a Bag-like potential 



V(r) 



V , if r>R; 
0, otherwise. 



(4) 



We also consider the one dimensional case for the nicety 
of the analytical result, 



V(z) 



if z < -f or z > f; 



V , 

0, otherwise 



(5) 
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The Vo — > oo limit gives, of course, the MIT Bag 
model 0, if one does not consider the Bag energy. 
We will consider models of vector+scalar confinement 



FIG. 1: 3-D Scalar Central Confinement 



FIG. 2: One Dimensional Scalar Square Well Confinement 



next, using in one case a linear spatial potential and on 
the other case a harmonic one, 



V(r) = ±(l + 1 °)(V + \r), 



V(r) = -(l + 7 °)Cr 2 . 



(6) 



(7) 



Finally we shall consider a purely scalar harmonic po- 
tential, 



V(r) = Cr 2 



(8) 



In the following sections it will be assumed that the 
current quark mass of 5 — 10 MeV can be neglected com- 
pared to the quark eigenenergy of several hundred MeV. 



II. SCALAR SQUARE- WELL POTENTIAL 

For any given potential V , from the Dirac equation we 
have that 



therefore 



— m)xf> = Vip, 



to) = (V). 



(9) 



(10) 



In the three dimensional case, for the central potential 
V(r) presented in (Q), the solutions of the Dirac equation 
for the ground state, with to = 0, in two regions I. r < 
R, and II. r > R (Fig. |l|) have the following form 



v I Jo(Er) , (s) 
ipi(r) = Ni . „ . )x w , 
\ia-rji{Er); 



ip II (r) = N 



ii 



h (ik r) 



(11) 



C is \ (12) 



where ko = \/V 2 — E 2 , and jo, ji are the spherical 
Bessel functions, and Hq , are the spherical Hankel 



functions of the first kind. The ground state energy can 
be found from the energy eigenvalue equation 



ji (ER) = j (ER) 



l + k R 



(V + E)R 



(13) 



while for Vq — > oo the eigenvalue equation is ji(ER) — 
jo(ER), as is familiar from the MIT Bag model [0 0. 

However, we know there are long range forces between 
baryons. If one wants to accommodate long range forces 
in this type of model, then one has to allow quarks to 
penetrate the walls of the potential well with some finite 
probability. Therefore the height of the potential, Vo, 
should be finite. A reasonable choice for Vo and R can 
be obtained by fitting the model parameters to obtain 
reasonable values, for example, for the mean square of 
charge radius of the nucleon (r 2 ) and for the axial vector 
coupling constant gA- We get a good fit by choosing 
R = 1.12 fm and V = 3 GeV for which we find (r 2 ) = 
0.64 fm 2 and gA = 1.15, as compared to experimental 
values of 0.76 fm 2 and 1.27 respectively Jl|. Solving (|l|) 
for this choice of parameters for the ground state energy 
of a quark we find E = 348 MeV. 



Using solutions given in (Till) and (12), we find 



- to) = 21 MeV. 



(14) 



Exploration of the integrals appearing in (V(r)), shows 
that (tf— to) — > as 1/Vb, when Vo — > oo. 

It may be of some pedagogic value to give the equiva- 
lent result for the ID case (Fig. ||). The wave function for 
\z\ < a/2 is just the free solution of the Dirac equation, 
and the solutions for \z\ > a/2 are obtained from the free 
solution by the substitution E — > E — Vq- We obtain 



: - to) = 2Vr 



ip ip dz 



E 



1 + a^V 2 - £ 2 



(15) 



One can note immediately that when the height of the 
potential Vo oo then (ji — to) — > 0, unless a — > 0. For 
the choice of parameters made above, we obtain 



: - to) = 14 MeV, 



(16) 



where for the ground state energy E we have used a value 
of 260 MeV, from the energy eigenvalue equation. 
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III. SCALAR + VECTOR LINEAR 
CONFINEMENT 

Let us consider now the confinement problem of a spin 
1/2 particle in a confining potential of the form 



V(r) = i(l + 7 °)(Vo + Ar). 



(17) 



This linear potential model for quark confinement was 
used in |D| to calculate several properties of low-lying 
baryons. In |pl| the authors assumed nonzero quark 
masses. The straightforward modification of the wave 
functions for the case of vanishing current quark masses 
yields the following solution for the lowest energy eigen- 
state of the Dirac equation for the potential (|l7|), 



*(r) = N[ ^ )x {s) 
\a-p/E $(r) ,X 



(18) 



$(r) 



h --Ai{Kr + ai )), (19) 



A-KAi l2 {al) r 



where K = {XE) 1 ^. The energy eigenvalue E and the 
normalization constant ./V are given in (po|) 



E = V -^, N 2 = 3E . 

K ' 4E - Vq 



(20) 



In [|13f| an analytic expression was obtained for the mean 
square charge radii of the barions and in Jl4[ Ferreira 
obtained an analytic expression for the magnetic moment 
of the proton. We modified those expressions for the 
zero current quark mass case and used them together 
with the energy eigenvalue equation < J20| ) to fit our model 
parameters Vq and A. We choose Vq = —626 MeV and 
A = 0.98 GeV/fm to fit (r 2 ) exactly and give the closest 
to the data value of /i p , obtaining 



E = 420MeV, 



= 0.76 fm 2 and fi p = 2.44 n.m. 



For the above mentioned values of the model parame- 
ters we find that 



27 MeV. 



(22) 



IV. SCALAR + VECTOR HARMONIC 
CONFINEMENT 

Consider now potential of the form 

V(r) = ^(l + 7 °)Cr 2 . 



(23) 



The solution of Dirac equation with this potential is given 
in p5| . They write the lowest energy state Dirac spinor 

as 



i>(r) = 



1 



4tt yr • fg(r)/r 



(24) 



where is a Pauli spinor, with the normalization 
J ip^ipd 3 r — J °°{f 2 + g 2 )dr = 1. Then the upper and 
lower components of the solution are 



fir) 

9{r) 



N 



r 

Tq 

N_ /r_ 

V3 \r 



-r 2 /2r 



(25) 



N = /sToW^), r 2 E 2 Q 



3, C 



:E* 



where Eq is the ground state eigenenergy and tq is a state 
dependent scale parameter. 

Now we can calculate the matrix element of interest, 



(tf-m) = I </ii(l +7 )Cr 2 ^ d 3 r 



C\ 1/3 E 



f(r) 2 Cr 2 dr 



(26) 

So, we can see that, in case of scalar+vector confinement 
of equal strengths, {ft — m) is determined only by spin 
independent part of the Dirac spinor and is equal to one 
third of the ground state energy. In jl(| it was also shown 
that for three massless quarks in their lowest Is orbit, 
with energy eigenvalues Eo for each quark, the center- 
of-mass energy obtained with the potential ( p3| ) is just 
£0, hence the nucleon mass in this model is: Afjy = 2i?o 
(instead of Afjy = 3Eq , as in non-relativistic and non- 
recoil models). Therefore, E — 540 MeV and, 



180 MeV. 



(27) 



V. PURE SCALAR HARMONIC POTENTIAL 

Tegen fl5|| has considered scalar+vector harmonic con- 
finement in calculating the weak neutron decay constant 
9a/ gv and found too small a value for gA/9Vi com- 
pared to experiment. In [ fl6| and [ jl7| , a pure scalar 
harmonic potential V(r) = Cr 2 was studied numerically, 
and yielded more satisfactory results for and for the 
RMS charge radius. We find that 



= C 



r 2 (f(r) 2 -g(r) 2 ) dr, (28) 



where f(r) and g(r) are defined as in eqn. (f24|). 

We have fitted the numerical solution presented graph- 
ically in [yj with C = 830 MeV/fm 2 to calculate our 
integral of interest (pq). The fitted wave functions are 
presented in Fig. |^, and as a benchmark for evaluation of 
the quality of the fit, we have calculated (r 2 ) and gA and 
obtained values 0.61 fm and 1.26 respectively, as com- 
pared to (r 2 ) = 0.64 fm and gA = 1-26 found in fjq| . 

Thus we obtained, without any additional tuning, the 
following result: 



(jf - m) = 160MeV. 



(29) 
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FIG. 3: Fit to the numerical solution of Dirac equation for 
pure scalar harmonic confinement. 

VI. SUMMARY 

In this paper we have calculated, for the ground state 
of the quark in a nucleon, the matrix element of the oper- 
ator (y-m), using variety of confinement potential mod- 
els, under the assumption that the constituent quarks 
obey the Dirac equation. The motivation has been to 
solidify the estimates of the non-commutativity parame- 
ter of canonical (Lorentz violating) noncomutative QCD, 
where some of leading order Lorentz violating effects are 
proportional to factors of (ji — to). 

Interestingly, we found the following results, 

!21 MeV, for 3-D scalar central pot. 
27 MeV, for scalar+vector linear pot. 
180 MeV, for scalar+vector harm. pot. 
160 MeV, for pure scalar harmonic pot. 

(30) 



We note that in the case of scalar central confinement 
as considered in section II, (rf — to) vanishes as 1/Vo 
when Vq — > oo, but it is different from zero in general. 
We note also that the value of — to) obtained for the 
scalar+vector linear confinement model is close to that 
obtained for scalar 3D potential well. 

We have also shown that in case of scalar+vector har- 
monic confinement of equal strengths, (j/— to) is deter- 
mined only by spin independent part of the Dirac spinor 
and is equal to one third of the ground state energy. 

For pure scalar harmonic confinement of the form 
V(r) — Cr 2 , — to) was obtained using a fit to the nu- 
merical solution of the Dirac equation presented graph- 
ically in [ p^[ , and appears to have a value pretty close 
to that obtained for the scalar+vector harmonic confine- 
ment model. 

Results obtained in this paper are within an order of 
magnitude agreement with the estimate made by Carl- 
son et al. The results obtained in [Q were used there 
to constrain the noncommutativity parameter in Lorentz 
violating noncommutative field theories. The constraints 
were very strong, and are still quite severe even if weak- 
ened by an order of magnitude. These results may be 
taken as a motivation to look for space-time noncomuta- 
tivity in Lorentz-covariant ways [JT^Sj, |l| 
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